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This paper presents a comparison between two approaches to sensor calibration. According 
to one approach, called explicit, an estimator compares the sensor readings to reference 
readings, and uses the difference between the two to estimate the calibration parameters. 
According- to the other approach, called implicit, the sensor error is integrated to form a 
different entity, which is then compared with a reference quantity of this entity, and the 
calibration parameters are inferred from the difference. In particular this paper presents the 
comparison between these approaches when applied to in-flight spacecraft gyro calibration. 
Reference spacecraft rate is needed for gyro calibration when using the explicit approach; 
however, such reference rates are not readily available for in-flight calibration. Therefore the 
calibration parameter-estimator is expanded to include the estimation of that reference rate, 
which is based on attitnde measurements in the form of attitude-quaternion. A comparison 
between the two approaches is made using simulated data. It is concluded that the performances 
of the two approaches are basically comparable. Sensitivity tests indicate that the explicit filter 
results are essentially insensitive to variations in given spacecraft dynamics model parameters. 


I. Introduction 

I N-flight sensor calibration is a crucial element in spacecraft (SC) preparation for a successful mission and much work 
has been done on this subject, the latest being the work reported in Ref. 1. During the first stage of the calibration 
process the sensor errors sources, which are called calibration parameters, are estimated. During the second stage, 
errors caused by these calibration parameters are removed from the sensor readings. 

There are two approaches to sensor calibration; namely, explicit and implicit approaches. In the explicit approach an 
estimator compares the sensor readings to reference readings, and it uses the difference between the two to estimate the 
calibration parameters. On the other hand, in the implicit approach the sensor error is integrated to form a different 
entity, which is then compared with a reference quantity of this entity, and the calibration parameters are inferred from 
the difference. If, for example, the sensors are accelerometers, then using the explicit approach the estimator compares 
the acceleration readings to reference acceleration. If die implicit approach is applied, then the accelerometer readings 
are properly integrated to generate velocity and the latter is compared with the reference velocity. The estimator uses 
this difference to estimate the calibration parameters. Similarly, if the sensors are gyroscopes, then in the explicit 
approach the estimator compares the measured rate with the reference one. In contrast, according to the implicit 
approach, the gyro readings are properly integrated to form attitude , which is then compared by the estimator to the 
reference attitude, and the estimator uses this difference to estimate the calibration parameters. As described, in both 
cases the estimator estimates the calibration parameters. Obviously, reference rate is required in the explicit approach to 
gyro calibration, and reference attitude is needed when the implicit approach is used. (It should be noted that by attitude 
we mean attitude in die wider sense; thus, two vector measurements constitute attitude reference in the wider sense 
because two vector measurements in the reference and body coordinates determine the attitude.) 

When gyro calibration is performed on the ground it is easy to measure the reference rates to which the gyro outputs 
are compared. In-flight, however, the reference rate is not readily available; therefore, the rate has to be estimated. This 
adds an additional complication to in-flight gyro calibration when the explicit approach is used. 

An example of the explicit approach can be found in Ref. 2. An improvement is presented in Ref. 3 where the 
calibration of the AQUA satellite gyroscopes was treated. Since the calibration was done in-flight, a reference angular 
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rate was not readily available. Therefore measured attitude was used to estimate the reference angular rate vector. 
Subsequently an estimator used the difference between the estimated reference rate and the gyro-measured rate to 
estimate the calibration parameters. Actually, the two estimators were combined into one Kalman filter. Because the 
rate estimation was based on Euler’s equation that describes the SC angular dynamics, the rate estimation was 
nonlinear. Classically, an Extended Kalman filter is used to estimate the state vector of a system in which either the 
dynamics or the measurement model or both are nonlinear. However, it was found in several cases 2 3 ' 5 that when Euler's 
equations describing SC angular dynamics are used in the dynamics model the use of the Pseudo-Linear Kalman Filter 
(PSELIKA) 6 produced very good estimates of the SC angular velocity. 

A simple block diagram that explains the explicit approach to calibration is shown in Fig. 1. A more detailed 
description is presented in Fig. 2. The implicit approach is normally used for calibrating the inertial measuring unit in 
inertial navigation systems 7 and is being used in spacecraft sensor calibration too (see e.g Ref. 8, 9). The methodology 
that governs the implicit approach when applied to gyro calibration is presented in a simple block diagram in Fig. 3 and 
a more elaborate description is shown in Fig. 4. 



Fig. 1 : Explicit Estimation of the Calibration Parameters. 
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Fig. 2: A Detailed Description of the Explicit Estimation of the Calibration Parameters. 
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Fig. 3: Implicit Estimation of the Calibration Parameters. 
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Fig. 4: A Detailed. Description of the Implicit Estimation of the Calibration Parameters. 


The purpose of this work is to compare between the explicit and the implicit approaches using simulated as well 
as real gyro data. In the next section we derive the gyro error model due to the calibration parameters. In Section in 
we develop the model employed by the explicit gyro-calibration estimator, and in the Section IV we present the 
model used by the implicit gyro - calibration estimator. Then, in section V, we show how the estimated calibration 
parameters are used to compensate the gyro readings. In Section VI we discuss some implementation issues, and 
show simulation results in Section VII. In the last section we draw our conclusions from this work. 

EL. The Gyro Error Model 

The gyro errors that are considered in this work are: misalignment, scale factor error, and bias (constant drift 
rate). The gyro error model is a linear model, which associates small error sources to the gyro outputs. Due to the 
linearity of the model we can model the contribution of each error source independently and then sum up all the 
contributions into one linear model. We start the description of the error model by deriving the expression for the 
gyro misalignments. 

II.l Misalignment Model 

The assumed direction of the sensitive axis of gyro x, which is one of the three considered gyro axes, is 
presented in Fig. 5. The body coordinate axes are also presented and are denoted by X, Y, and Z. The orientation of 
this gyro is expressed by a vector of unit length in the direction of the gyro sensitive (input) axis. It can be seen that 
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Fig. 5: Misalignment definition of the x-gyro sensitive axis. 


the resulting rate error due to misalignment of all three gyros is: 
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where HI - is the projection of the i-gyro sensitive axis on the j body-axis that is perpendicular to the i body-axis. 
Since the sensitive axis is described by a unit vector and due to the proximity of the gyro sensitive axis to the 
respective body-axis, m- is a misalignme nt angle, assumed to be small, expressed in radians. The elements CD W 
are the angular velocity components measured by the gyros. Let, 
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Then Eq. (1) can be rewritten as. 
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Ato" = Q"m 


(4) 


II.2 Scale Factor Error Model 

As mentioned, additional error sources that cause the difference between the correct value of the actual 
rates and their meas ure ments are the scale factor errors. The error model for the scale factor errors is simply 

co k 

z x 

to k (5) 

(i> k 

z z J 

where the superscript k denotes die fact that this error is caused by gyro scale factor errors, assumed to be small, and 
kj is the scale factor error of gyro i, i = X, y, Z . Eq. (5) can be written as follows 
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113 Bias Model 


The bias error model is quite simple and is given by 

Ata b = I 3 b 


( 8 ) 


where I 3 is the third order identity matrix and 


b 

x 



(9) 


where x, y, z are the corresponding gyro axes. 

n.4 The Augmented Gyro Error Model 

The total gyro error is the sum of all the error discussed before; namely, misalignment, scale factor and bias 
errors; that is 
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A© = A©“ + A© 1 + A(o b 


Using Eqs. (4), (7.c), and (8) we get 


(10 a) 


A© = f2 m m + £2 k k + I 3 b 

The last equation can be written in the following form 


(10 b) 


A© = [q“ Q k I 3 ] 


Define H(o r ) ((O r is the real angular rate vector) as follows 



H(©)=[fT Q k I 3 ] 

(lO.d) 

Also let 


x T = [m T k T b T J 

(lO.e) 

then Eq. (lO.c) can be written as 


A© = H(© r )x 

(10.f) 
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ID. The Explicit Gyro Calibration Algorithm 

Attitude information and rates measured by the gyros are available in-flight but, unlike ground calibration, the 
reference rates needed for explicit calibration are not readily available. Therefore, as mentioned earlier, in this case 
we estimate the angular rate vector while estimating the calibration parameters. We use the attitude information to 
estimate the angular rate. The attitude information can be supplied in various ways; namely, it can be in die form of 
raw vector measurements or it can be given in an already processed form as attitude quaternion, for example. The 
estimation of the angular rate vector hinges on the following Euler’s equation that describes the SC angular 
dynamics 

© = r'[(l©+h)x]©+r 1 (T-h) (11. a) 


where I is the SC inertia tensor, [(ItO + h)x] is the cross product matrix of the vector (Ico + h) where a cross 
product of a general vector a is defined as follows 
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h is the angular momentum of the momentum wheels, h = dh / dt and T is the external torque operating on the 
SC. 

As mentioned before, the reference rate is estimated from the measured attitude. This necessitates the inclusion 
of the attitude kinematics in the estimator dynamics equation. Suppose that the attitude is measured by autonomous 
star trackers (AST) 10 that yield the measured quaternion. The quaternion kinematics equation is expressed by 11 



(12.a) 
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Equation (12.a) can be also written as 


(12.b) 


where 



(12.c) 


(12-d) 


Because X is a constant vector, it obeys the following differential equation 

i = 0 (13) 

We can combine Eqs. (11), (12.c) and (13) into one dynamics equation, add white noise to the angular dynamics 
equation to account for model uncertainty, add a small quantity of white noise to the calibration parameter dynamics 
to better model them in case they are not really constant but rather vary slowly*, and add white noise to the 
quaternion kinematics equation to account for modeling inaccuracies. As a result we obtain the following dynamics 
equation 
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where is the zero-mean white noise vector added to the angular dynamics, W x is the zero-mean white noise 
vector added to the calibration parameters, and W q is the white noise added to the quaternion dynamics. Because the 
dynamics matrix is a function of <0 as well as q, the dynamics equation presented in Eq. (14) is nonlinear, 
subsequently a nonlinear estimator is needed for estimating the state vector. The most appropriate estimator for the 
nonlinearity structure of this dynamics model is the PSELIKA algorithm 3 . 

As explained before, we have two measurements to consider; namely, the gyro rate-measurements, <D m , and the 
quaternion measurement, q m . The former consists of the true rate, the error due to the calibration parameters, X , 
and some high frequency noise which is modeled as a zero-mean white noise, V ra . Consequently we have 


a m = to + Aco + v o (15.a) 

Using Eq. (lO.f), the last equation can be written in the following matrix form. 


* It is good practice to add some white noise to states in the filter model even when they are known to be constant (see Ref. 12). 
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(15-b) 
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The quaternion measurement is modeled by a combination of the true quaternion and a zero-mean white 


measurement noise vector, V q . 
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As explained in the introduction, we combine the estimator that estimates the calibration parameters and that which 
estimates the rate vector into one estimator; therefore, we have to properly combine the two measurement equations. 
The merger of Eqs. (15.b) and (15.c) yields 
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Equation (14) and Eq. (15.d) constitute respectively the dynamics and the measurement equation for die augmented 
estimator used in the explicit approach to gyro calibration. 

IV. The Implicit Gyro Calibration Algorithm 

As in the previous case, our goal now is to estimate X , and, following the logic behind the implicit estimation 
approach, for that we need to know how X influences the attitude estimation. The true quaternion; that is, the 
quaternion that describes the hue SC attitude, solves the differential equation (12.a) where Q is a function of the 
true angular rate vector, to , which we do not know. We know and therefore we would rather use, the measured rate 

vector, to . According to Eq. (15.a), to = CO + AtO + V m therefore 
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and where 6Q and 5Q V are functions of CO and \ a respectively, and are in the format of Eq. (16.c). Using Eq. 
(16.b), Eq. (12.a) can be written as: 
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which can also be written as 
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(17 b) 
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where Q is given in Eq. (12.d). Using Eq. (lO.f) in Eq. (17.b) yields 
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Augmenting the last equation with Eq. (13) to which, as before, we add W x , yields 
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Eq. (17.d) is the dynamics model used by the estimator when we apply die implicit approach to gyro calibration. The 
measurement is the attitude measurement (and not the gyro measurements), which in our case is simply tj m . From 
Eq. (15.c) it is obvious that here 


q„=[i 



(18) 


Equation (17.d) and Eq. (18) constitute, respectively, the dynamics and the measurement equation for the estimator 
used in the implicit approach to gyro calibration. 


V. Compensation 

To complete the calibration process we need to perform its second stage; namely, compensation of the gyro 
readings using the estimated calibration parameters. This stage is performed when either die explicit or the implicit 
approaches are used. It can be performed periodically during the estimation process or at its end. Compensation 
during the estimation process is preferable because then the models used by the filters, which are functions of the 
hue angular velocity, become more accurate as the estimation proceeds. From Eq. (lO.c) we obtain 
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A calibrated gyro measurement, CD , is computed as follows 
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Using Eq. (19.a) in the last equation yields the final gyro compensation formula 
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(19.d) 


VI. Implementation Considerations 

The dynamics matrix in Eq. (14) is a function of © and q which are not available. Initially, however, we can 
evaluate the matrix using q m rather than q but when the attitude converges we can use the estimated quaternion, 
q , which is expected to be closer to q . Similarly, initially we can use CO m , or the compensated gyro outputs (when 
compensation is done during the estimation process) rather than © r (which is not available to us) to evaluate the 
dynamics matrix of Eq. (14), the measurement matrix, H , of Eq. (15.d), and the dynamics matrix in Eq. (17.d). 
When the estimator converges in the explicit approach, it is better to switch to the estimated rate denoted by © , and 
when the estimate of X converges in the implicit approach, it is better to use the calibrated rate, © , rather than © m 
because the calibrated gyro readings are normally expected to be closer to © r more than © m js. 


VII. Test Results 

The algorithms were tested using simulated flight data. An angular rate profile was selected which assured 
complete observability of the estimated calibration parameters (see Fig. 6). It was assumed that no external torques 
operated on the SC. Thus, T was set to zero. In order to apply the correct wheel momentum that generated the 
desired ©profile, Eq. (ll.a) was transformed into 

h = [-©x]h +[(I©x)]© - 1© + T (20) 


and was solved for h where © and © were, of course, known functions of time. The correct attitude quaternion 
was generated using the prescribed © . The associated AST readings were simulated by adding the noise vector, 
V q , to the measured quaternion (see Eq. (18)). That zero-mean noise was generated using a Gaussian random 

number generator with a standard deviation of 10 5 . The following calibration parameters were used in the 
simulation. 


T 

m = 

[ 0.1 

- 0.1 0.1 

- 0.1 

0.1 - 


~o.of 



~+ 0 .f 

k = 

0.01 

( 21 . b ) 

b = 

- 0.1 


0.01 



+ 0.1 


deg/sec 


(21. a) 


(21 -c) 


The simulated gyro readings were generated using © , the calibration parameters, and a random noise vector, V ffi 

(see Eq. (17.c). The latter vector was generated using a zero mean Gaussian random number generator with a 
standard deviation of 0.8 deg/hr. The nominal inertia tensor was: 
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100 0 0 
1 = 1 0 200 0 
0 0 300 i 


kg-m 2 


(21. d) 


Each sensor provided data at a 1 Hz rate. The attitude profile started with an inertial period of 200 seconds (see 
Fig. 6) during which practically all of the gyro biases were estimated. This inertial period was followed by a 0.1 
deg/sec maneuver about the x-axis that lasted 200 seconds, followed by another 200 seconds of inertial period Then 
at 600 seconds a 0.1 deg/sec maneuver about the y-axis was initiated which, too, lasted 200 seconds and was 
followed by 200 seconds of inertial period. At 1000 seconds a 0.1 deg/sec maneuver commenced, this time about file 
z-axis, lasting also 200 seconds. When the latter maneuver ended at 1200 seconds, die SC stayed inertial until the 
end of the run at 1400 seconds. The two algorithms were applied to the AST readings and the algorithms produced 
estimates of the calibration parameters. In order to compare the performance of the two algori thms, figures of merit 
were computed for the misalignment, scale factor, and bias estimates based on the final values of the calibration 
estimates. The figures of merit were defined as follows: 
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where: 


£ mi = i* misalignment final estimation error, 
ntj = 1 th misalignment. 

8 k ; = i ft scale factor final estimation error, 
kj = i* scale factor. 

£ b j = i* bias final estimation error, 
bj = i* bias. 

VII. 1 Results with the Explicit Filter 

Figure 6 presents the true angular rates during the simulation (in black) and the estimates (in red). We note that 
good rate estimates are imperative for obtaining good calibration parameter estimates when using the explicit 
approach. The true and estimated values of the gyro misalignments are presented in figure 7. The scale factors are 
compared in figure 8, and figure 9 presents the same for the gyro bias values. 


VII.2 Results with the Implicit Filter 

The true and estimated values of the gyro misalignments are presented in figure 10. The scale factors are 
compared in figure 1 1, and figure 12 presents the same for the gyro bias values. 
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Fig. 9: True (black line) and explicit estimate (red line) of gyro biases. 
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Fig. 12: True (black line) and implicit estimate (red line) gyro biases. 


V1I.2 Comparison of the Figures of Merit 

The three pairs of figures of merit, one for the misalignment, one for the scale factor and one for the bias 
estimation errors are presented in Table I. It is seen that while the explicit algorithm es tima tes the misalignments and 
scale factors better than the implicit algorithm, the latter does better in estimating the biases. However the changes 
are minor and depend heavily on our ability to tune each filter. In the scenario presented, the performance of the two 
approaches is comparable. 


Table I: Figures of Merit 



Explicit 

Implicit 

J m (for misalignment) 

0.14987098 

0.19447241 

J k (for scale factor) 

0.07091744 

0.01163662 

J b (for bias) 

0.00033426 

0.00028284 


VII.3 Sensitivity Tests 

Unlike the implicit filter, the explicit filter uses a model of the SC dynamics in terms of Euler’s equation. In the 
preceding simulations of the explicit filter we used in the filter-model the SC parameters that were also used in the 
simulation of the angular rate profile (we designate this case as the nominal case). It is important to examine the 
sensitivity of the filter results to uncertainty in these parameters. Table I presents the comparison between the 
nominal explicit nominal filter and the implicit filter. In Table II we present a comparison between the results 
obtained when using the nominal explicit filter and those obtained when the inertia matrix in the filter is 90% of its 
nominal value. In table Id we present a comparison between the nominal case and the results obtained using the 
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explicit filter where 1^ in the filter is 10% higher than its nominal value, I is 10% lower, and 1^ is higher than 

its nominal value. Finally, in Table HI we present the comparison when the components of the filter h are 10% 
lower than their nominal values. It is obvious that the filter is not sensitive to parameter uncertainty. 

Table II: Figures of Merit of the Nominal Case and the Case of 10% Reduction in the Inertia 



Nominal 

Mismatched 

J m (for misalignment) 

0.14987098 

0.14968341 

J k (for scale factor) 

0.07091744 

0.07080843 

J b (for bias) 

0.0 0033426 

0.00033426 


Table III: Figures of Merit of the Nominal Case and the Case of Alternate Inertia Variation 



Nominal 

Mismatched 

J m (for misalignment) 

0.14987098 

0.15031518 

J k (for scale factor) 

0.07091744 

0.07087014 

J b (for bias) 

0.00033426 

0.00033426 


Table IV: Figures of Merit of the Nominal Case and the Case of 10% Reduction inh . 



Nominal 

Mismatched 

J m (for misalignment) 

0.14987098 

0.14989170 

J k (for scale factor) 

0.07091744 

0.07091747 

J b (for bias) 

0.00033426 

0.00033426 


VHL Conclusions 

This paper presented a comparison between two approaches to die problem of SC sensor calibration. In 
particular, the problem of gyro misalignment, scale factor, and bias estimation was considered. The first approach, 
die explicit approach, utilizes attitude measurements to obtain rate estimates, which are compared to the gyro- 
measured rates in order to estimate the calibration parameters. The rate estimation and the calibration parameter 
estimation are performed by one pseudo-linear Kalman filter. The second approach, the implicit approach, uses the 
un-calibrated gyro outputs to compute an attitude that is compared with the measured attitude. The difference, which 
is a function of the calibration parameters, is also used in a pseudo-linear filter. In this paper the error models used 
by the filters of the two approaches were developed and excessive simulation runs were made for different values of 
calibration parameters. It was found that the performances of the two approaches were comparable. However, the 
results are sensitive to filter tuning. Sensitivity tests were carried out on the explicit filter dynamics model 
uncertainty. It was found that the explicit filter results were practically insensitive to inertia and wheel momentum 
uncertainties. In the present simulations quaternion observations were used, which were based on star tracker 
measurements. It is recommended that sun sensor and magnetometer measurements, which provide less accurate 
attitude determination, also be tested. 
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